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We present the analysis of a nonlinear effect of parametric 
oscillatory instability in power recycled LIGO interferometer 
with the Fabry-Perot (FP) cavities in the arms. The basis for 
this effect is the excitation of the additional (Stokes) optical 
mode with frequency lu\ and the mirror elastic mode with fre- 
quency u m , when the optical energy stored in the main FP 
cavity mode with frequency too exceeds the certain threshold 
and the frequencies are related as uo — coi + u m . The pres- 
ence of anti-Stokes modes (with frequency u\ a — luq + cu m ) 
can depress parametric instability. However, it is very likely 
that the anti-Stokes modes will not compensate the paramet- 
ric instability completely. 



I. INTRODUCTION 

The full scale terrestrial interferometric gravitational 
wave antennae are in process of assembling and tuning at 
present. One of these antennae (LIGO-I project) sensi- 
tivity expressed in terms of the metric perturbation am- 
plitude is projected to achieve the level of h ~ 1 x 10~ 21 
[1,2]. After the improvement of the isolation from noises 
in test masses (the mirrors of the 4 km long optical FP 
cavities) and after increasing the optical power circulat- 
ing in the resonator up to W ~ 830 kW the sensitivity 
is expected to reach the value of h ~ 1 x 10~ 22 [3]. This 
value of W corresponds to the energy £ — 22 J stored in 
the FP resonator. 

In previous paper [4] we have described the possibly 
existing effect of pure nonlinear dynamical origin which 
may cause substantial decrease of the antennae sensitiv- 
ity or even the antenna disfunction. The essence of this 
effect is classical parametric oscillatory instability in the 
FP cavity which modes are coupled with mechanical de- 
gree of freedom of the mirror. This effect appears above 
the certain threshold of the optical energy £q when the 
difference luq — uii between the frequency uq of the main 
optical mode (which stores So) and the frequency uj\ of 
the idle (Stokes) mode is close to the frequency oj m of 
the mirror mechanical degree of freedom. The coupling 
between these three modes appears due to the pondero- 
motive pressure of the light photons in main and Stokes 
modes and due to the parametric action of mechanical os- 
cillation on the optical modes. Above the critical value of 
energy £q (dimensionless parameter TZo > 1, see below) 
the amplitude of mechanical oscillation will rise expo- 
nentially as well as the optical power in the idle (Stokes) 
optical mode. 

In the article [4] we have used a simplified model of this 
effect in which lumped model of mechanical oscillator had 



been used and optical modes with gaussian distribution 
over the cross section had been taken into account. Under 
these assumptions the parametric oscillatory instability 
will appear if 
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Here m is the value of the order of mirror mass, L is 
the distance between the FP cavity mirrors, Au>i — luq — 
u>i — Lj m , Si and Qi = u)\/25\ are the relaxation rate 
and quality factor of the Stokes mode correspondingly, 
S m and Q m — oj m /28 m are relaxation rate and quality 
factor of the mechanical oscillator. 

Recently E. D'Ambrosio and W. Kells [5] have reported 
that if in the same one dimensional model the anti-Stokes 
mode (with frequency uj\ a = ujq + uj rn ) is taken into ac- 
count then the effect of parametric instability will be sub- 
stantially dumped or even excluded. In this article we 
present the analysis based on the model which takes into 
account several important details of the antenna. This 
analysis shows that parametric oscillatory instability still 
may exist. 

In section II we present the analysis of one dimensional 
optical model of antenna in which the so called power 
recycled mirror is taken into account. In section III im- 
portant key elements of 3-dimensional approach are used 
to prove that it is very likely that the anti-Stokes modes 
will not compensate the oscillatory instability. 



II. THE ROLE OF POWER RECYCLING 
MIRROR IN THE ANTENNA 

The design of laser interferometer gravitational wave 
antenna apart from the two main optical FP cavities 
also includes the so called power recycling mirror (PRM) 
which allows to increase the value of £q using the same 
laser input power (see fig. 1). 
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FIG. 1. The scheme of power recycled LIGO interferometer 
with the Fabry-Perot (FP) cavities in the arms. The end 
mirror of the FP I cavity is a part of mechanical oscillator. 
T bs and T pr are transmittances of the beam splitter and power 
recycling mirror correspondingly, li is the distance between 
the input mirror of the FP I cavity and the beam splitter, 
l pr is the distance between the power recycling mirror and 
the beam splitter, Fi and Ei are the complex amplitudes of 
travelling optical waves in different parts of interferometer. 

The results of calculations presented below are based 
on the following simplifying assumptions: 

• Only one mirror (in the FP I cavity in fig. 1) is 
movable, it is a part of mechanical oscillator that 
is considered as lumped one with single mechani- 
cal degree of freedom (eigenfrequency w m , quality 
factor Q m = iv m /25 m and mass m which is of the 
order of the total mirror mass). All other mirrors 
are assumed to be fixed. 

• All mirrors have no optical losses. Both end mirrors 
of FP cavities have the ideal reflectivity. The input 
mirrors of the FP cavities are identical and have the 
finite transmittances T\ = T 2 = T = 2n L / '(XoQ op t) 
(Aq is the optical wavelength, Q op t is the quality 
factor, L is the distance between the mirrors). The 
distances L\ = Li = L are equal. 

• Only the main (pumped) mode with frequency u>o 
and relaxation rate 5o = u>o/2Qq, and Stokes mode 
with lo\ and b\ = LO\j2Q\ correspondingly are taken 
into account (Q and Qi are the quality factors). 
It is assumed that ujo — uJi — uJ m - 

• Only the main mode is pumped by laser and the 
value of stored energy £ is constant (approxima- 
tion of constant field). 

• In this particular model we do not take into account 
the possible influence of the anti-Stokes mode. 



It is possible to calculate at what level of energy £ 
the Stokes mode and mechanical oscillator become un- 
stable. The origin of this instability can be described 
qualitatively in the following way: small mechanical os- 
cillations with resonance frequency ui m modulate the dis- 
tance L that causes the excitation of optical fields with 
frequencies loq ± u> m . Therefore, the Stokes mode ampli- 
tude will rise linearly in time if time interval is shorter 
than relaxation time. The presence of two optical fields 
with frequencies w an d Wi will produce the component of 
ponderomotive force (which is proportional to the square 
of fields sum) at difference frequency — w\. Thus this 
force will increase the initially small amplitude of me- 
chanical oscillations. 

For analysis of parametric instability we have to use 
two equations for the Stokes mode and mechanical oscil- 
lator, and find the conditions when this "feedback" pre- 
vails the damping which exists due to the finite values 
of relaxation in mechanical resonator and in the FP cav- 
ity. Note that in resonance w — + <^m the effect of 
parametric instability for power recycled LIGO interfer- 
ometer is larger than for the separate FP cavity because 
the Stokes wave (at frequency uj\) emitted from the FP 
cavity throughout its input mirror is not lost irreversible 
but returns back due to power recycling mirror, therefore, 
its interaction is prolonged. 

We can write down the field components Eq", E\ n of 
the main and Stokes modes inside the FP I cavity corre- 
spondingly and the displacement x of mechanical oscilla- 
tor in rotating wave approximation as: 



E™ = A [D Q e-^ ot + D*e Wot ], 
E[ n = A^Dxe^ 1 * + Dle i0J % 
x = Xe~ lUmt +X*e lUm \ 



(3) 



where Do and D\ are the slowly changing complex am- 
plitudes of the main and Stokes modes correspondingly, 
and X is the slowly changing complex amplitude of me- 
chanical displacement. Normalizing constants Ao , Ai are 
chosen so that energies So, l stored in each mode (of the 
FP I cavity) are equal to £ , l = w o, llA), i | 2 /2. Then one 
can obtain the equations for slowly changing amplitudes 
(see details in Appendix A): 
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{dt + 6 m ) x = 



m(jj m L 



(5) 



Remind that we assume D to be constant. The addi- 
tional relaxation rate 6 pr describes the relaxation of os- 
cillations in the FP cavity with power recycling mirror. 

One can find the solutions of (4, 5) in the following 
form 
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£>!(*) = £>Ie 
X(i) =Ie 



Ait 



Ai = A + iAuJi 



and write down the characteristic equation: 
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The parametric oscillatory instability will appear if one 
of the characteristic equation roots real part is positive. 
Analysis of this equation with assumption 



S m *C 5 pr <C 5i, 



gives the condition of parametric instability (see details 
in Appendix B): 



TZ 



2 + 



Si 
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> 1. 



(8) 



In ultimate resonance case when Au>i <C S pr we obtain 
the condition of resonance parametric instability: 



TZ x 1 + 



Si 



25 



pr 



> l. 



(9) 



It means that the parametric instability takes place at 
energy £ which is smaller than one for the separate FP 
cavity by the factor of <~ ^T^- 

It is important that factors TZo and have large nu- 
merical values, using the parameters for probe mass fab- 
ricated from high quality [8] fused silica that are planned 
to use in LIGO-II we obtain (see [6] and Appendix D): 



TZo ^ 6100, 



A 

25 



31 . 



(10) 



pr 



It means that if Au>i <C 5 pr the maximal energy stored in 
the FP cavity can not exceed the value of about 6100 x 
31 ~ 1.9 x 10 5 times smaller than one planned for LIGO- 
II (!). 

Note that the estimate of TZo differs from the estimate 
presented in [4] because here we use parameters more 
close to ones planned for LIGO-II. In particular: (a) mass 
is four times greater than in [4] ; (b ) mechanical frequency 
is assumed to be two times smaller; (c) loss angle (and 
mechanical relaxation rate) is ~ 4 times smaller in accor- 
dance with results of Ageev and Penn [8]; (d) relaxation 
rate Si is ~ 6 times smaller. 

The probability of the ultimate resonance (Awi <C Spr) 
is extremely low because the value of 5^ 



pr 



1.5 s" 



1S 



rather small. In more realistic case when Aui\ ^> 5iS pr 
we obtain the " partial resonance" condition of parametric 
instability from (8): 



TZo 



(>+¥) 



> l 



(ii) 



It means that in the case when Si » Awi 3> y/SiS pr the 
maximal energy stored in the FP cavity can not exceed 
the value of about 3050 times smaller than one planned 
for LIGO-II (!). Note that condition (11) differs from the 
parametric instability condition in the separate FP cavity 
[4] by the factor of 2 in denominator. 

The considered model is the simplest one and more 
detailed 3-dimensional model of interferometer has to be 
analyzed. In particular, there are reasons to hope that 
the danger of parametric instability may be smaller in 
real interferometer: 

• Even in resonance the Stokes and elastic modes may 
not spatially suit to each other (small overlapping 
factor) . 

• The possible presence of the anti-Stokes mode may 
partially or completely depress the parametric in- 
stability. 

In the next section we consider both these factors. 



III. CONSIDERATION ON 
THREE-DIMENSIONAL ANALYSIS 

It is possible to generalize the above simplified model 
for the arbitrary elastic mode in the mirror. It has been 
shown [4] that in this case the constant TZo m condition 
(8) should be multiplied by the overlapping factor Ai. 

In general case when both Stokes and anti-Stokes 
modes have to be taken into account the characteristic 
equation may be presented in the following form: 



(X + 5 rn ) = TZ SiS m x 



Ai 



Si 



Xi+Si [ 2(Ai 
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2(Aia + Spr a ) 

V (J fo(r±)fi(r±)u z dr ± ) 2 
f\fo\ 2 d?±J\fi\ 2 dr ± J\u\*dV 

V (J Mr±).fia(r ± )u z dr ± ) 2 
f\f \ 2 dr ± f\fi a \idr±f\u\*dV> 

A + iAwi, Au>i = u>o — u>i — Lo-m , 

X + iAlVia, AlVi a — Ld la - LU - U) r , 



The equation (12) is the generalization of (6). Hereinafter 
subscript i corresponds to the Stokes mode and subscript 
ia to the anti-Stokes mode. Ai, Ai a are the overlapping 
factors for the Stokes and anti-Stokes modes correspond- 
ingly, /o, fi and fi a are the functions of the optical fields 
distribution in the main, Stokes and anti-Stokes optical 
modes correspondingly over the mirror surface. Vector u 
is the spatial vector of displacements in the elastic mode, 
u z is the component of u normal to the mirror surface, 
/ df± corresponds to the integration over the mirror sur- 
face, and / dV over the mirror volume V. 
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Analyzing this characteristic equation one can see that 
the presence of the anti-Stokcs mode can considerably 
depress or even exclude parametric instability. For ex- 
ample, let the main, Stokes and anti-Stokes modes be 
equidistant and belong to the main frequency sequence 
wi = n(K - l)c/L, uj = ttKc/L, u la = n(K + l)c/L 
(K is an integer). In this case Alu\ — Ac<Ji a , the main, 
Stokes and anti-Stokes modes have the same gaussian dis- 
tribution over the cross section and hence the same over- 
lapping factors: Ai — A la . It means that the second 
term in the right part of (12) is larger than first term, 
the positive damping introduced into elastic mode by the 
anti-Stokes mode is greater than negative damping due 
to the Stokes mode, hence the parametric instability is 
impossible. This case has been analyzed in details in [5]. 

However, it is worth noting that this situation is possi- 
ble only for the small part of the total number of op- 
tical modes (see fig. 2). Indeed, resonance conditions 
uo ~ ivi + uo m can be fulfilled with a relatively high prob- 
ability for many of the optical Stokes and mirror elas- 
tic modes combinations. If wc assume the main optical 
mode to be gaussian with the waist radius of the caustic 
wo (the optical field amplitude distribution in the mid- 
dle between the mirrors is <~ e~ r ^ w °)), then the Stokes 
and anti-Stokes modes are described by generalized La- 
guerre functions (Gauss-Laguerre beams) and the set of 
frequency distances between the main and Stokes (anti- 
Stokes) modes is determined by three integer numbers: 



uj - wi 



2(2A + M) L\ 
K arctan 



27TWq 



L \ 7T 

~ (2.4 if -0.66 A- 0.33 M) x 10 5 s _1 , (13) 
uj ^ (2.4 K a + 0.66 N a + 0.33 M a ) x 10 5 s _1 . 



where Ao is the wavelength, K = ± 1, ±2 ... is the lon- 
gitudinal index, N = 0, 1,2..., and M — 0, 1,2... are 
the radial and angular indices, other numerical parame- 
ters are given in Appendix D. We see that full depression 
of parametric instability takes place only if the Stokes 
and anti-Stokes modes belong to the main sequence i.e. 
M = N = 0, when Awi = Awi a and the Stokes and anti- 
Stokes modes have equal spatial gaussian distribution. 
Such modes obviously present the small part of the total 
optical modes number. Indeed for K = 1, N = M = 
we have from (13) ojq — o->i — 2.4 x 10 5 sec -1 . However, 
our numerical calculations show (see below) that the low- 
est elastic mode has the frequency of about nine times 
smaller: w m lowest — 0.28 x 10 5 sec -1 and only within the 
range between 0.28 x 10 5 sec -1 and 1.6 x 10 5 sec -1 there 
are more than 50 (!) elastic modes and each of them 
has to be carefully considered as possible candidate for 
parametric instability. 

For the case when the Stokes and anti-Stokes modes do 
not belong to the main sequence (non zero numbers N 
and M) the frequencies of the suitable Stokes and anti- 
Stokcs modes are not equidistant from the main mode 
(i.e. Awi ^ Awia) and have different spatial distribu- 
tions (i.e. Ai 7^ Ai a ). Illustration of this is given in fig. 2. 
For the shown Stokes mode (left to the main mode) there 



is no suitable anti-Stokes mode (it should be located right 
to the main one). In this case one can use the approxi- 
mate condition for the parametric instability (see details 
of approximations in Appendix C): 



Au 



(14) 



Aw 



> 1 



Here the second term in the right part describes the in- 
fluence of the anti-Stokcs mode. 



Main sequence (M = N = 0) 



^_ Main mode 



LAJ 



U) 



FIG. 2. Schematic structure of optical (Laguer-Gauss) 
modes in the FP cavity. The modes of the main frequencies 
sequence are shown by higher peaks. It is shown that Stokes 
mode with frequiency wi may not have suitable anti-Stokes 
mode (it is denoted by question-mark). 

The above consideration for lossless mirrors can be gen- 
eralized for mirrors with losses. Analyzing the more im- 
portant case when only the FP cavity mirrors have losses 
with loss coefficient R < T for each of them we can easily 
show all previous formulas to be valid taking the following 
substitutions into account: 

2R 

6 1 =>5 1r = 5 1 (1 + ti), ??= — <!, (15) 



<>, -r S prR = ( -jp +T] ) /), . 



(16) 



We see that in the case r] <§C Tp r /4 all formulas does not 
change, in the case of 1 S> rj 3> T pr /A the condition of the 
ultimate resonance changes but the condition of partial 
resonance remains unchanged. For LIGO-II the losses are 
negligible (it is planned that R ~ 5 ppm, r/~2x 10~ 3 <C 

Tpr/A). 

The account of losses is important when we consider 
optical modes with high indices M, N because their 
diffractional losses increases for higher indices M, N. We 
have numerically calculated the equivalent loss coefficient 
Rm,n describing the diffractional losses on each mirror 
and have found that 



2i?M. n T, 
Vdif = < 



~ 0.015, if 2N + M < 6, 
r] dii <l, if 2A + M<9, 
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We see that there is a wide range of indices M, N which 
must be taken into account when analyzing parametric 
instability. Note that even for the case rjdif > 1 the para- 
metric instability may also be possible due to the large 
numeric value of factor IZo- 

Using the Femlab program package we have numer- 
ically calculated the first 50 elastic modes frequencies 
and spatial distributions in the cylindric mirror fabri- 
cated from fused silica (parameters are listed in Appendix 
D). These frequencies lied within the range between 
28000 sec" 1 and 164000 sec" 1 . We have estimated the 
errors Aw m of these frequencies comparing the data ob- 
tained with different total number of nodes (N n ~ 1135 
and N n ~ 3777): Aw ra ~ 500 ■ ■ ■ 4000 sec" 1 (the error 
increases for elastic modes at higher frequencies). Such 
error is unacceptably large because in order to determine 
the ultimate resonance in power recycled interferometer 
the error have to be Aw m <C S pr ~ 1.5 sec -1 or in the 
case of "partial resonance" Aw m <C #i — 94 sec -1 . 

It is worth noting that the frequency density of elastic 
modes rapidly increases at higher frequencies. In partic- 
ular, the mean distance Au m d between the elastic modes 
frequencies can be estimated as follows: 



Aw 



md 



m lowest 



and it is equal to Aw m j ~ 100 sec -1 ~ <5i even at u> m = 
3.8 x 10 5 sec -1 . It means that practically for each elastic 
mode with frequency higher than 3.8 x 10 5 sec -1 there 
exists the Stokes mode with small detuning: Aui < 5± 
("partial resonance"). However, the same speculations 
can also be referred to the anti-Stokes modes and hence 
the accurate calculation of overlapping factors Ai, Ai a 
is required for each of the elastic modes. 

Even in the case of parametric resonance the overlap- 
ping factor Ai may be zero (for example, elastic mode 
and the Stokes mode can have different dependence on 
azimuth angle). However, it is important to take into 
account that only the elastic mode attached to the mir- 
ror axis in contrast to the optical mode which can be 
shifted from the mirror axis due to non-perfect optical 
alignment. Hence, the overlapping factor must depend on 
distance Z between the center of mirror and the center of 
the main optical mode distribution over the mirror sur- 
face. It means that Ai may be zero for Z — but nonzero 
for Z 0. Therefore, the numerical analysis of the mode 
structure should evidently include the case when Z^O. 
Note that there is a proposal to use special shift Z of the 
laser beam of about several centimeters from the mirror 
axis in order to decrease thermal suspension noise [9] . 

Due to the necessity to decrease the level of thermoe- 
lastic and thermorefractive noises [10,12,11,7] the size of 
the light spot on the mirror surface is likely to be substan- 
tially larger and the light density distribution in the spot 
is not likely to be gaussian (the "mexican hat" modes 
[7]) to evade substantial diffractional losses. The optical 
modes which are complementary to such a "mexican hat" 
main mode have the frequencies more close to the main 
mode than given by equation (13), and the probability 



to be entrapped into the parametric instability is higher. 
Thus the estimates presented above for gaussian optical 
modes may be regarded only as the first approximation 
in which the use of analytical calculations is still possible. 



IV. CONCLUSION 

Summing up the above calculations and considerations 
we have to conclude that the effect of parametric oscil- 
latory instability is a potential danger for the gravita- 
tional wave antennae with powerful laser pumping. From 
our point of view, to estimate correctly this danger it 
is necessary to implement the following subprograms of 
researches: 

1 . To calculate numerically the values of eigenfrequen- 
cies uj m for lower elastic modes with relative errors 
of at least 10~ 4 •••10~ 3 . It is not an easy task 
because the error in standard schemes of finite ele- 
ments calculations rises as square of u> m . It is nec- 
essary to keep in mind that these calculations will 
play a role of introductory ones because it is likely 
that in LIGO-II non-gaussian mode distribution of 
light will be used (so called "mexican hat" mode) 
and, correspondingly, it will be necessary to cal- 
culate numerically all spectrum of "mexican hat" 
modes. 

2. At the same time the numerical analysis can not 
solve the problem completely because the fused sil- 
ica pins and suspension fibers will be attached to 
the mirror. This attachment will change the elastic 
modes frequency values (and may be also the qual- 
ity factor and distribution). For example, assuming 
the pin mass of about Am ~ 80 g one can estimate 
that frequency shift may be about 



A < A- 

mpm ~ 1x105 sec^ 1 



Am , „„„ -i 
x - — ~ ±100 sec" 1 , 
2m 



i.e. about the value of Si (!). In addition, the 
unknown Young modulus and fused silica density 
inhomogeneity (we estimate it may be of about 
10~ 3 • • • 10~ 2 ) will additionally limit the numeri- 
cal analysis accuracy. Thus we have to conclude 
that the direct measurements of eigenfrcquency val- 
ues, distribution and quality factors for several hun- 
dreds of elastic modes for each mirror of FP cavity 
are inevitably necessary. 

3. When more "dangerous" candidates of elastic and 
Stokes modes will be known their undesirable influ- 
ence can be possibly decreased. Perhaps, it can be 
done by accurate small change of mirror shape or 
by introducing low noise damping [13]. 

4. The last stage of this program should be presented 
by the direct tests of the optical field behavior with 
smooth increase of the input optical power: it will 
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be possible to register the appearance of the pho- 
tons at the Stokes modes and the rise of the Q m in 
the corresponding clastic mode while the power W 
in the main optical mode is below the critical value. 

We think that the parametric oscillatory instability effect 
can be overcome in the laser gravitational antennae after 
these detailed investigations. 
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5i — ifl 
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APPENDIX A: POWER RECYCLED 
INTERFEROMETER 

In this Appendix we deduce the equations (4, 5). 

1. Internal FP Cavities 

For the fourier components of complex slowly changing 
amplitudes in general case we have obvious expressions 
(see notations in fig. 1): 



Ft n (Q) 



2i5 1 F 1 (Q) 



£i(fi) = -Fi(fi) 



Ti (Si — iAuj — itl) 
Si + iAu + iCl 



81,2 



4Li 2' 



Si — iAu — ifl 



(Al) 



(A2) 



One can obtain the similar formulas for the second FP 
cavity: 



Fi n (n) 



E 2 (Q) = -F 2 (Q) 



2iS 2 F 2 (ty 

/%(s 2 -iny 
s 2 + in 



5 2 — id' 



(A3) 
(A4) 



2. The Mean and Small Amplitudes 

Let us introduce the mean amplitude and small ampli- 
tude (denoted by small letters). For example, for ampli- 
tude in the FP I cavity it means: F[ n = F[ n + f{ n . 

Let us also assume that 5 2 = Si, and now we keep in 
mind that 



Auj 



UJlX 



Then the mean amplitudes are written as follows: 
2i 



2i 



Ti, Ei — —Ti, 
F2, £2 = —J~2 ■ 



(A5) 

(A6) 
(A7) 



3. Beam Splitter 



Assuming that Tb s = 1/2 one can write the following: 

F ie-^ = -L ( t F 3 +F 4 ), F 2 = ^= (F 3 + zF 4 ) , 

Ei = -L ( lE2 + Ei e** 1 ) , E 3 = ±= (E 2 + iE 1 e** 1 ) 

Here (f>i — kli is the wave Ei phase shift due to length 
path li between the FP cavity and the beam splitter. We 
assume that <j>i — n/2 and analogous phase shift 4> 2 (for 
the wave E 2 ) is equal to zero. Then one can obtain for 
the mean amplitudes the following: 



-t— T 3 

rr- ?3 -T 3 

£4 = 0, £3 = —T 3 , 

Small amplitudes are equal to: 

, -h + ih f h + ifi 
h — 7= — , J2 - 



V2 



V2 



(All) 

(A12) 
(A13) 



(A14) 



r r , i^Ti j-31. 1 
e4 = h I H 2 — ' 3 = ~-' 3 ~ ' 



4. Power Recycling Mirror 

We have the following expressions: 



= i^%-F 5 + y/l-TprEae*", 
E 5 = i^% r E 3 e 1 ^ + ^/l-T pr F 5 , 

(uj + ALO pr + Q) 



(A15) 
(A16) 

(A17) 



Using (A13) and assuming that the PR cavity is in reso- 
nance: exp(ioj l pr I c) = i (i.e. Uol pr /c = ir/2 + 2nn, n is 
an integer) one can obtain: 



-2T b F 2F 5 

J~3 — /Tf^ l ^3 



IT 

' 1 rt 



hA 3 



T3T1 



(Ais) 

(A19) 



B SOLUTION OF THE CHARACTERISTIC EQUATION (8 ) 



A 3 ~ i-r + r - 



1 - 



i(Aw pr + n) 



(A20) 



cT 



pr 



up r 



The value of 8 pr » | /\uo pr + fi| and hence one can rewrite 
(A20) as follows: 



A 3 i - r 

x Tp r Sl 



Tn 



o- 



<5 pr - if? 
<5i — iSl 



"pr - 4 • 

And finally, for / 3 one can obtain the following: 



h 



Si — iCt 



2(S pr — in) 
Using (A16)we will have: 
£5 = 



e 5 



T pr fo^o + 



+ /5 



(A21) 
(A22) 

(A23) 

(A24) 
(A25) 



However, such substitution will be incorrect, because the 
equation for the fourier transform of complex amplitude 
has not contain time-dependent terms. The correct form 
of this equation is the following: 



/f n (fi) {6 pr -in) 



(A31) 



iD™X*{n - Awi) 



- in + ■ 



Now one can obtain time domain equivalent of (A31) 
applying inverse fourier transform: 



(a t + *i)(St + <f pr )/i n (t) = 

= iDp e - iAuJlt — x 
£1 



iAwi 



(A32) 
X*(t) = 



= iD, 



Zi" 



<9( + <5pr + 



6l 



X*(t)e 



The last equation is the complex conjugate to (4) where: 

/r -> Di. 

7. Equation for Elastic Oscillations 



5. The Small Variational Amplitude in the FP I 
Cavity 

In order to calculate the small amplitude fi one can 
rewrite equation (Al) using (A14) and (A19) assuming 
that small amplitudes and f$ describing vacuum fluc- 
tuations of the input waves are zero: 



{Tl n + /f) x (S 1 - iAuj - in) = 

_ iv^fi (-JT 3 - f 3 + if A ) 



1 

2iT x i^2T 3 _ 2i^2T h 

\/T\Fp T 



/i n (Ji -ifi) =«Awf; n 
iV2 6x 



£1 - in 
2(<v - in) 



(A26) 

(A27) 
(A28) 
(A29) 



It can be rewritten using (A27) as: 



fl n (5i - in) (S P r - in) = iJ^Aco 



^n -|- 



For elastic displacement we have the following equa- 
tion: 



d t x + 2S m d t x + u m x = 
_2S(£o+£i) 2 ^ 



m 



Air 



ire V 1 

1)4 ), 



(A33) 



(A34) 



n\* rin i{ujQ—uj 



where S is the cross section of light beam, c is the light 
speed. Introducing the slow amplitudes for displacement 
x: 

x(t) = X{t) e"^™ 4 + X*{t) e 1 "™' , 
one can obtain the following: 
iS 



dtX + 8 m X 



■in f rin\* „ — iAu>it 



2ircmtjj. 
Aui = u>o - wi - uj rn . 



(A35) 
(A36) 



(A30) 



This equation coincides with (5) after the substitutions, 
listed below: 



fir. 
Jl 



Di, To -» Do . 



6. Time Domain 



Now it seems that one can make in (A30) the following 
substitution: 

Tf -> Doe-*'*'-" 1 * , Aw^^ X *e iuJmt 



APPENDIX B: SOLUTION OF THE 
CHARACTERISTIC EQUATION (8) 

In this Appendix we obtain the instability condition 
(8) from the characteristic equation (6). 



C SOLUTION OF THE CHARACTERISTIC EQUATION (14) 



8 



Let us write down the the solution of equation (6) as 
sum of real and imaginary part: 

A = a + ib . 

The condition of instability is a > 0. Thus substituting 
A = ib into (6) one can find two equations (introducing 
notations A and B): 



U Si5 m x 



(b + A Wl ) 2 + 6 2 



(Bl) 



We write down the the solution of this equation as a 
sum of real and imaginary parts anew: 

A = a + ib 

and assume a = 0. So substituting A = ib into (12) 
and extracting real and image parts we can write two 
equations: 



S m = H 5i5 m x 



*iAi 



{b + A^f + 5 2 



(CI) 



1 + 



diSp r -(b + Awi) 2 



2[(b + Auj 1 r + Sn 



1 + 



SiSpr - {b + Aujx) 2 



2[(6 + A Wl ) 2 + e] 



b = -K 5iS m x 
x 1 + 



(b + Auj^ + S 2 

Sl(Sl + 5 pr ) 



2[(b + A^) 2 + S 2 



(B2) 



— K 1<V> X 



Slavic 



Wi (6 + A^ la ) 2 + <5 2 a 



A la 



SlgSpra - (o + Acj la ) 2 

2[(b + Au la ) 2 + <5 2 ] 



J m (fe + A^QC 
5iB 



(B3) 



Using the last equation one can formally express b as: 
5-m.C / / LC 



-Awi 



5-i.B 



1 



5iB 



C _ 25 2 pr + 5 1 (h + 5 pr ) + 2(b+ Aujxf 



B 2S%\ + + (b + Awi) 2 

It is obvious from the last expression that 



(B4) 



2 < — < -r^- 



B 



jpr 



for any value of (6 + Awi ) . Hence one can conclude from 
(B4) that |6| <C Awi (remind that o~ m <C o~i) and equation 
(Bl) can be simplified as follows: 



S m = K 5iS m x 



(A^f + S 2 ) 
S pr (Si + 2<5 pr ) + Alu 2 
+ 



(B5) 



Now rewriting this equation one can easily obtain the 
condition of parametric instability (8). 



APPENDIX C: SOLUTION OF THE 
CHARACTERISTIC EQUATION (14) 

In this Appendix using some approximation we de- 
duce the instability condition (14) from the characteristic 
equation (12). 



AiBi-AiaBla, 



i „ , , Ai(& + A Wl ) 

6 = -7£ <Mm X : - — X 



(b + A^Y + 5 2 

Sl(8\ + Spr) 



2[(b+Aco 1 ) 2 +S 2 ] 



+ 



Ci 



i T> £ £ Wla Ala (6 + Awi a ) 

+ft <Mm x — - — x 

wi (6 + Auj la ) 2 + Sf a 

5la(5l a + 5 pra ) 



2[{b + Auj la ) 2 +5 2 a ] 



= -A X C, 



Ci 

6 + Awi 



+ A\ a C\ a 



b + Aui a 



(C2) 



(C3) 



(C4) 



Ol Ola 

Using notations (C2) this equation can be rewritten as: 
6 = — t — j5 — x (6 + Awi) + -j— - — x (6 + Awi a ) , 

Ol±>l Ola-Dla 

and one can formally express 6 as: 

I —At,}-, <?mlC ' 1 4- A(,)i Sml aCl a 1 



= 



(C5) 



Using definitions (CI, C3) it is easy to prove that 

2 < < -i-L or "^ ml < ( ^ m1 ^' 1 ^mi 
-Bi <5p r <5i S\B\ Sp r 

C\a _ 0~l o 2o~ m i a 5ml aC\ 

2 < — — < or - - < — — < 



B\a &pr a 



8 a l Si a B la d pra 
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Now assuming that 

6 -f- « 1, ^ « 1 , (C6) 

Op r Op r a 

we can conclude that 

\b\ < Awi, Awi„. (C7) 

Then the parametric instability condition (14) can be eas- 
ily obtained from (CI). 

It is worth noting that inequalities (C6) are the key 
assumption for deduction of the PI condition (14). These 
inequalities obviously correspond to the case when values 
S ml = A\B\ and 8 m ia — A\ a B\ a are not very close to 
each other. In order to understand how close they can be 
let us assume that 

Smla = 8 m i(l - e), e<l, 

and try to estimate thcminimal value of e. From (C6) one 
can obtain 8 m \ ~ S m ia — S m /e. Hence the inequalities 
(C6) are equivalent to: 

4^ « 1 • (C8) 

Using parameters of LIGO-II (see Appendix D) we have 
estimates 6 m ~ 1 x 10" 3 sec -1 , d pr ~ 1 sec -1 . Hence one 
can conclude that for values of e > 10 -2 the inequalities 
(C6) fulfill. 
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APPENDIX D: NUMERICAL PARAMETERS 

For interferometer we use parameters planned for 
LIGO-II. More details see in [6]. 



(Dl) 



Here W is the power circulating inside FP cavity (£q = 
2LW y We assume that cylindric mirror (with radius R, 
height H and mass m) is fabricated from fused silica with 
angle of structural losses <^ = 1.2xl0 -8 [8]: 





= 2 x 10 15 sec" 1 , 


w 


= 5.5 cm, 


T 


= 5 x 10 -3 , 


L 


= 4 x 10 5 < 


T 


= 6 x 10 -2 , 




~ 10 m, 


Si 


~ 94 sec -1 , 




~ 1.5 sec - 


w 


= 830 kW, 


So 


~ 2.2 x 10 : 



R = 19.4 cm, H = 15.4 cm, 

m = 40 kg, p = 2.2 g/cm 3 , 

E = 7.2 x 10 11 cr = 0.17, 

4> = 1.2 x 10 -8 , 5 m = uo m cj)/2. 



(D2) 



Here E is Young's modulus, a is Poison ratio. It is useful 
to calculate factor TZo for these parameters and elastic 
mode frequency iv m = 1 x 10 5 sec -1 : 

TZ — ^l^m _ £0^1 ^ g 

2mL 2 uj^ l SiS m mL 2 uo 2 n (t>5i 



